clear;
clc;
A=vander([1 2 3 4 5 6 7]),
a=poly(A);
B=polyvalm(a,A);
e=norm(B)

disp('----------------------')
a1=poly1(A);
B1=polyvalm(a1,A);
norm(B1)